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Abstract 

When a planar shock wave interacts with a random pattern of pre- 
shock density non-uniformities, it generates an anisotropic turbulent 
velocity/vorticity field. This turbulence plays an important role at 
the early stages of the mixing process in the compressed fluid. This 
situation emerges naturally in shock interaction with weakly inhomo- 
geneous deuterium- wicked foam targets in Inertial Confinement Fusion 
(ICF) and with density clumps/clouds in astrophysics. We present an 
exact small-amplitude linear theory describing such interaction. It is 
based on the exact theory of time and space evolution of the per- 
turbed quantities behind a corrugated shock front for a single-mode 
pre-shock non-uniformity. Appropriate mode averaging in 2D results 
in closed analytical expressions for the turbulent kinetic energy, de- 
gree of anisotropy of velocity and vorticity fields in the shocked fluid, 
shock amplification of the density non-uniformity, and sonic energy 
flux radiated downstream. These explicit formulas are further sim- 
plified in the important asymptotic limits of weak/strong shocks and 
highly compressible fluids. A comparison with the related problem 
of a shock interacting with a pre-shock isotropic vorticity field is also 
presented. 

1 Introduction 

The interaction of shock waves with flow inhomogeneities attracted the at- 
tention of scientists working in different fields, ranging from Aerodynam- 
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ics [SI El El El [SI HOI HSl HH], Shock tube research [JJ [21 E], Laser Fusion 
[121 HSl HZl [221 [23 ESI I2H1 EDI Ell E21 E31 [Ml ESI ES] and Astrophysics 
[371 138| 139] . The physics of the interaction has been described theoretically 
0HHm[3[Tni[IIl[T2]and experimentally [231 M, 122] . In the last decades, 
the advent of super computers provided a sophisticated and accurate tool 
with which to follow the details of that interaction [121 ITB I H^ l H51 fTE \ 1201 
I2T| 122] . However, the development of analytical works is an important and 
complementary way to guide in the design of numerical simulations and ex- 
periments. This work is a natural continuation of a previous article [T7J, in 
which the interaction of a planar shock with a pre-shock turbulent vorticity 
perturbation field was considered, with the help of an analytical theory. The 
shock/density field interaction has been studied in the recent past by differ- 
ent researches with the help of numerical simulations. Rotman [12] provided 
results of that study assuming 2D perturbations upstream of the shock and 
Mahesh et al [13] studied a full 3D spectrum of steady state density per- 
turbations. With the help of 2D Large Eddy Simulations (LES), Rotman 
described the interaction of a shock wave with random vorticity and density 
fields separately. For both cases, he noticed that shock compression waves 
greatly amplifies the upstream turbulence and reduces the turbulent length 
scales [12]. On the other hand, Mahesh et al have used Direct Numerical Sim- 
ulations (DNS) to study the interaction of a combined spectrum of pre-shock 
vorticity and density inhomogeneities. Besides, the problem considered here 
is of importance in the field of Inertial Confinement Fusion (ICF), where a se- 
quence of well tuned shocks must be launched inside the thermonuclear foam 
targets. The use of foams has been proposed over 20 years ago to improve the 
stability of the implosion and give an effective smoothing of the laser energy 
deposition process [221 EDI EH E2J E31 EH E51 EE]- However, given the non 
uniform character of the background density through which the first shock 
travels inside the target, it is clear that the shock front will dynamically re- 
act against those perturbations, generating additional density, vorticity and 
acoustic fluctuations in the compressed fluid, much in the same way as hap- 
pens in the shock/vorticity interaction studied in [T7J and references therein. 
The anisotropic turbulent spectrum left by the shock front can act as a trigger 
for further hydrodynamic instabilities and mixing. In this work we present a 
completely analytical model to study the outcome of the interaction of a pla- 
nar shock front with a non uniform pre-shock 2D density profile, continuing 
the task started in [T7J. The strategy is simple: an arbitrary pre-shock spec- 
trum may be decomposed in its Fourier modes and the shock interaction with 
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Figure 1: A planar shock moving with velocity Dx in the laboratory reference 
frame hits a 2D turbulent density field located in the half-space x' > 0. 

every single mode is studied separately and superposition is later performed 
to see the effect on the whole spectrum. Each shock/single mode interaction 
is assumed to be linear, that is, the initial density fluctuation is assumed to 
be much lower than the mean background density, a condition satisfied in 
ICF targets. Statistical averages are obtained by integrating over the angles 
that define the orientation of the pre-shock perturbation wavenumber vec- 
tor in space, similarly as has been done in Refs.pjU [17] and the works cited 
there. We will restrict our results to isotropic pre-shock spectra. This work 
is structured as follows: in the following Section 2, the 2D single mode theory 
is developed, and we show the exact shock dynamics with its corresponding 
asymptotic expressions. In Section 3 the 2D isotropic random density field is 
considered. We obtain exact expressions for the kinetic energy and vorticity 
generation and density amplification, and acoustic energy flux emitted by 
the shock. 

2 Interaction of a planar shock with a single 
mode 2D density field 

2.1 Wave Equation and Boundary Conditions 

A planar shock is incident at t — at the surface x' = in the laboratory 
frame. The fluid is an ideal gas with adiabatic exponent 7 and the pertur- 
bations consist of a weakly inhomogeneous density field. The shock comes 
from the left (x' = —00) and travels with velocity Dx', as measured in the 
laboratory reference system. In the uniform half-space x' < 0. the density 
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and pressure ahead the shock are pi, pi respectively, and p 2 , P2 are the val- 
ues behind it. The velocity of the compressed fluid is Ux', also measured 
in the laboratory frame reference. The upstream sound speed is c\, and the 
downstream value is c 2 . The shock Mach number with respect to the up- 
stream gas is Mi = D/c\ > 1 and the shock Mach number with respect to 
the compressed fluid is M 2 = (D — U)/c 2 < 1. Before the shock arrives to 
the interface x' — 0, the relationship between the quantities at both sides the 
shock front are: 

p, _D_ (7 + l)Ml m 
p 2 D-U ( 7 -l)M 1 2 + 2' [) 



Mo = 



C'2 



(7 - I) Ml + 2 
c 2 ~ \] 2 7 M 2 - 7 + 1 ' 

P2 _ 2jM[ -7 + I 
Pi 7 + 1 

v /(2 7 M 1 2 -7 + l)[(7-l)M 1 2 + 2] 



(2) 
(3) 
(4) 



ci (7 + l)Mi 

In the right half-space, the perturbed density field is described by 8pi(x', y) = 
pi€k cos(k x x') cos(k y y) where x' and y are the longitudinal and transverse 
coordinates as measured in the laboratory system of reference. To remain 
within the limits of validi ty of th e linear theory, we assume < 1, where 
efc is a function of k = ^Jk^ + ky, as isotropy is assumed for the pre-shock 
perturbations. The longitudinal and transverse wave numbers are defined, 
respectively by: k x = 2ir/\ x and k y = 2ir/X y , with and X y characteristic 
lengths. Once the shock is in the half-space x' > 0, the density profile in front 
of it will induce density, pressure and velocity fluctuations downstream, and 
its shape will be distorted. We assume that all the perturbed quantities are 
much smaller than the background values. We denote the upstream values 
with the subscript 1 and the downstream values with the subscript 2. From 
now on, the linearized equations of motion will be written and solved in a 
system of reference that co-moves with the compressed fluid particles. We 
define the followings dimensionless perturbation functions, factoring out the 
small parameter e^: 

— = e k p(x,t) cos(k y y) , 
Pi 
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= e k p(x, t) cos(k y y) 

P2C 2 

e k v x (x,t) cos(k y y) 



c 2 

e k v y (x,t)sm(k y y) . (5) 



c 2 

In Eq.Q £ is the time, and x is the longitudinal coordinate as measured 
in the compressed fluid frame. The quantities v x and v y correspond to the 
longitudinal and transverse velocities respectively, and p and p represent 
the dimensionless density and pressure perturbations. We also define the 
dimensionless time r = k y c-it- The linearized equations of motion in the 
compressed fluid frame are: 

dp dv x 



dr d (k y x] 

dvr dp 



dr d (k y x) 

-- P, (6) 



8Vy 

Or 



representing the mass, x-momentum and y-momentum conservation equa- 
tions, respectively. Furthermore, the conservation of entropy holds if we 
assume adiabatic flow behind the shock, which is represented by: 

dp_ = dp_ 

The dynamics of the perturbed quantities in the whole compressed fluid is 
governed by the wave equation. Combining Eqs.fjS]) and ((Tj) we get the wave 
equation for the pressure fluctuations: 

dr 1 = d{k y x) 2 ~ P ' ( ' 

We assume that the shock front hits the interface x = at t = 0. At t — + , a 
shock is transmitted to the right into the perturbed half-space and a neutrally 
stable sound wave is reflected back inside the region x < 0, as shown in Figj2j 
As the shock wave travels in a non-uniform fluid, the shock will be distorted. 
We define the shock ripple ip s {y, t) as the deviation from planarity. The shock 
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ripple oscillates in time, generating pressure fluctuations that propagate with 
the local sound speed into the compressed fluid. The sound waves generated 
by the shock oscillation can be stable or evanescent waves, depending on the 
ratio k x /ky, the shock Mach number M\ and the gas compressibility 7 [T7] . 
We assume that no sound wave hits the shock surface from behind (isolated 
shock). At the surface x — 0, pressure and normal velocity are continuous on 
both sides of it. However, the distorted front generates vorticity and entropy 
perturbations and the neutral sound wave to the left does not. Hence, the 
rc-derivative of v y and density are generally discontinuous at x = 0. It is not 
difficult to get the following relationships just to the right of the left traveling 
sound wave: 



v x + P = 



. 



(9) 



The boundary conditions at the shock are obtained after linearizing Rankine- 
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Figure 2: (a) A planar shock wave travels with velocity D at the laboratory 
reference frame before enters into the perturbed density field . (b) The trans- 
mitted corrugated shock moves into the disturbed fluid at the compressed 
fluid frame for t > with velocity (D — U). Density, vorticity and acoustic 
fluctuations are generated behind it. A stable sound wave is reflected to the 
left traveling with velocity C2 in this frame. 



Hugoniot conditions and using the continuity of the tangential velocity. We 
write them here for the particular case of density /entropy pre-shock modula- 
tion at constant pressure and zero velocity perturbation ahead of the shock: 



(It 



Ml 
1M-M, 



7 + 1 - 
4M 2 P 

+ 1 _ 



M 2 R 5pi 



P- 



2 

M 2 (R 



Pi 

_l)5pi 
Pi 



(10) 
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v y = M 2 {R-l)i s . (12) 

where Eq. ffTUj) represents the mass equation, Eqs. ffTT]) and f[T2"l) correspond 
to the longitudinal and transverse momentum conservation respectively, and 
Eq.(!T3l) is the energy equation. Here, £, s €k = k y ip s is the dimensionless shock 
ripple amplitude. To get the perturbed quantities in the compressed fluid 
we solve the wave equation [Eq.©] with adequate boundary conditions. To 
this scope, we concentrate at the weak discontinuity (x = 0) and at the 
shock front [x = x s (t) = (D — U)t\. To finally solve the dynamics of the 
shock front, we need the initial value of the pressure perturbation behind the 
shock. Using Eqs.Q, ffTTT) and f[T3l) it is easy to get [T6] : 



MfM$(R - 1) 
Ps0 ~ 2M 2 M 2 + M 2 + V { } 

where p s o is the initial shock pressure perturbation at t = + . Besides, it is 
clear that the initial shock ripple amplitude is £ s o = 0, as the shock front is 
planar in shape when it arrives to x = 0. 

2.2 Pressure dynamics in the compressed fluid 

To solve the wave equation [Eq.©] inside the compressed fluid, we use the 
following coordinate transformation [U [161 El I2Z1 [28] : 

k y x = r sinh x, 

t = r cosh x- (15) 

Here, x = const represents a planar front defined by x = C2ttanhx. The 
shock front coordinate is given by: tanhx s = M 2 , and from Eq. (|T5|) we get: 

t = r s cosh Xs = —f= S (16) 
'I - M, 2 



The wave equation [Eq.flS])] is now rewritten as: 

,d 2 p dp 2 _ c 

h r h r p = - 

dr 2 dr dx 2 



2 d 2 p dp 2 _ d 2 p 
r l^ + r ^Z + r P= 1^- ( 17 ) 
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The boundary conditions at the shock front [Eqs. (110j) -f ll3j) ] can be recast as: 

ldp\ Ml + 1 dp, MliR-1) AUR-l), . u , 

~7T = - .)w2w 7, / = & o Co sin (CoM , 18 

rdxJxs 2M(M 2 dr s h _ M | 2 



^ = 7 / + 1 £ = cos (Cor s ) , (19) 

dr » 4MW1-M? 2,/l-M? 



where Co is a dimensionless frequency that characterizes the periodicity of 
the pre-shock density inhomogeneity. Its value is given by: 

k x M 2 R , s 

Co = ~~r~ • 20 

^ Jl - Ml 



We solve Eqs.( IT7|) -( ITU|) by using the Laplace transform. For any quantity 
<p(x, r ) we define its Laplace transform by: $(x, s) = <p(x, r ) exp(— sr)dr. 
After some algebra which we omit here, but that can be found in [T6l [T7] , we 
get an exact closed form expression for the Laplace transform of the shock 
front pressure fluctuations P s : 



p M 2 M 2 (R- l)s 

s ^ " 2MlM 2 s^Tl + (M? + 1) s 2 + Ml + 

2M 2 M 2 a e s 



2MlM 2 sVs T TT + (Ml + 1) s 2 + Ml\ (s 2 + Co) 
The coefficient a e , which is the same as a 2 o in [16] is given by: 



(21) 



a,. 



M 2 (R-1) f M 2 R 



1 -M 9 2 



Co 2 , (22) 



The evolution of the shock pressure perturbation as a function of the time 
is obtained by calculating the inverse Laplace transform of Eq. (l2~T]) . after an 
integration in the complex plane. We formally write [101 HTJ W2\ : 

rc+ioo _ 

Ps(Ts) = —.\ P s (s) exp(sr s )ds , (23) 

It is clear that the denominator of Eq. fl2ip contributes with stable asymptotic 
oscillations of frequency Co- After analyzing the residues at the poles s = 



8 



±i(o, we get the asymptotic expressions, similarity as in Ref . [T7] : 

e lr cos ( ( r s ) + e h sin ( ( r s ) , Co < 1 

(24) 

e s cos ( ( r s ) , , Co > 1 



where r s = J 1 — M| r, and the coefficients ei r , en and e s are the same as 
those obtained in [16], and formally equivalent to those shown in Eq.(47) in 

= 2M 1 2 M 2 [M 1 2 -(M 1 2 + 1)C 2 ] a e 
Qr 4MfM|Co 2 (1 - Co 2 ) + [Ml - {Ml + 1) Cof ' 



( : 'u 



4M?Mj( Jl - Co 2 a 



AMfM%Q (1 - Co 2 ) + [Ml - {Ml + 1) Cof ' 
2M 2 M 2 a e 

2MlM 2 ( Ja ~ 1 + (Mi 2 + 1) Co 2 - Ml ' 



(25) 



The shock front ripple oscillates with the dimensionless frequency Co within 
the domain of the variable r s . If Co > 1 sound waves fill the space behind it 
in the form of traveling fronts. At any position x, the sound waves will reach 
the asymptotic regime when the shock is far enough (x s ^> x). Because of the 
Doppler effect, the dimensionless frequency of the compressed fluid particles 
oscillations is Ci < Co- The value of Ci can be seen to be given by [161 fTT] : 



Co-M 2V /Co 2 -l 
Ci = , V '° 2 • (26) 



Besides, a longitudinal wave number k a ^ is associated to the stable fluid 
oscillations given by [THl HZ] : 



k = m ^q - vco 2 -i ^ (27) 

k v J I - Ml 



It is not difficult to see that the asymptotic pressure oscillation at any position 
x, is therefore given by: 

p(x, t) = e s cos (Cit - k a x c x) . (28) 
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For £o < 1, the sound waves emitted by the shock front are evanescent and 
decay exponentially away from it [H)J[T7]. From Eq. (l2"7|) it is easy to see that 

for 1 < Co < 1 / J 1 — Mf , the waves are emitted to the right, following the 

shock front, and for ( > 1 / J 1 — Mf the sound waves escape to the left, 
filling the whole compressed fluid. 



2.3 Rotational and irrotational perturbations down- 
stream 

As the shock ripple oscillates in time, not only pressure oscillations are gener- 
ated downstream, but also vorticity and velocity fluctuations. The vorticity 
is strictly generated at the shock front and remains frozen to the fluid ele- 
ments in the absence of viscosity. Let us define the 2D dimensionless gradient 
operator: 

^ f d d \ 



2D 



dk y x 



dk y y 



(29) 



It is clear that the vorticity is directed along the z-axis and is defined by: 



kyC 2 (V 



2D X V 



(30) 



with v = (v x ,v y ). The function uj can be found by using the continuity of 
tangential velocity at the shock front and that vorticity is conserved along 
the particles paths. Following [T7J |27] we get: 



kyxjl - Ml 



M 2 



+ fl 3 cos (Rk x x) 



where fl 2 arid ^3 are given by: 



sm(k y y) , (31) 



(M 1 2 -l) 2 V / 27M 1 2 -7 + l 
M 1 2 [(7-l)M 1 2 + 2] 3/2 

M 2 (R 2 - 1) 



n 



(32) 



(33) 



The first term, proportional to ^2 is generated by the shock front distortion 
(or shock curvature as named in [HI [19]). It is always present whenever the 
shock front gets corrugated as in any RMI like problem [T7]. The second 
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term originated from the interaction between the pre-shock density field and 
the zero order pressure jump across the shock surface, usually called the 
baroclinic term [181 US]- We omit a term Qi in Eq.( l3~Tj) that would appear 
for a pure pre-shock vorticity field ahead of the shock wave, a case that has 
been studied in detail in [17]. After combining Eqs.flS]) and (J7J), the velocity 
field downstream can be seen to satisfy the differential equation [T7J : 

^Id~v + V 2 d x (V 2D x ff) . (34) 



Or 2 

To solve it, we decompose the velocity into a rotational and acoustic compo- 
nent: 

v(x, y, t) = v rot (x, y) + v ac (x, y, t) . (35) 

The rotational contribution, which accounts for the vorticity downstream, is 
time independent in the compressed fluid frame, and it satisfies: 

Vl D V ro t = -V 2 D X (V 2 D X U rot ) , (36) 

and the acoustic part satisfies the homogeneous wave equation: 

= VL^c • (37) 

An exact solution can be sought for both contributions in the same way as 
has been done in [T7J for the pre-shock vorticity case. As we are interested 
in the asymptotic velocity field behind the shock, the asymptotic solution 
can be easily obtained, by matching the downstream asymptotic velocities 
with the asymptotic expressions at the moving shock wave. The procedure 
is similar to the calculations shown in [T7J and will not be repeated here. For 
the longitudinal rotational part we get: 

Q l rot cos (Rk x x - (fi rot ) cos (k y y), Co < 1 

v r f(x^X y ,y) 

Q s rot cos (Rk x x) cos (k y y), Co > 1 



and for the transverse component: 

R]rQ l rot sin (Rk x x - <p rot ) sin (k y y), Co < 1 

RbQrot sin (Rk x x) sin (kyy), Co > 1 



(38) 



f^(x»A,,y) 



(39) 
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where it is assumed that the point x is far enough form the interface x = 
and from the shock front [x <C (D — U)t]. The quantities Q l rot and Q s rot are: 



(°&rot — / , \ 2 ' ^rot — / , \ 2 V 4U / 

1 + ( R t) 1 + ( R k) 

and 

tan ^ = oJTsw ' (41) 

The acoustic velocity field is: 

v x c (x > \ y , y, t > 1) = Q ac cos (Cit - fc^z) cos (k y y) , (42) 



r > \ y , y, t > 1) = ^ sin (dr - k a x c x) sin (A;„j/) , (43) 



v Ci 



M>Co - \/Co - 1 
Q ac = V /° e s , (44) 



where: 



only valid for the short wavelength regime (Co > 1)- 

3 Interaction of a planar shock with a 2D ran- 
dom density field 

In this section, we study the interaction of an initially planar shock wave 
with a 2D isotropic random density field. We get averages of the downstream 
kinetic energy and vorticity generation. We also obtain the density amplifi- 
cation and the acoustic energy flux emitted by the shock. As we have seen in 
Section 2, the single mode profile is characterized by the vector k = (k x , k y ). 
Its components are: 

k x = kcosO , 
k y = -k sin 9 , (45) 

where we define k = \k\, and 9 is the incidence angle of the upstream pertur- 
bation wave vector with respect to the axis perpendicular to the shock front 
(2). The range of variation is < 9 < n. We consider an isotropic profile 
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and we assume that the wavenumber vector k is uniformly distributed along 
the unit semicircle, i.e. the probability of a particular orientation is dO/n. 
We can write the dimensionless frequency (o as a function of the incident 
angle as follows: 

M 2 R 1 

CO = ; „!.■■./!! ( 46 ) 



tan 6*1 



3.1 Turbulent kinetic energy generation 

In the previous section, we have obtained the asymptotic velocity profiles 
downstream for a single-mode pre-shock density field. Those profiles are 
useful to calculate statistical averages over a full spectrum of pre-shock non- 
uniformities. Let us introduce a dimensionless conversion coefficient K be- 
tween the pre-shock density non-uniformity and the post-shock turbulent 
Mach number (M turb ): 



MLb 



Kxe 2 k 



(47) 



where v 2 



v 2 

X 



v 2 , and for the 2D problem: 



K 



2D — — 



? r^/ 2 / 








IX Jo K 





de 



(48) 



The polar angle 9 can be expressed as a function of Co : 



d9 



MiJR(M 2 - 1) 



M 1 2 i?+(M 1 2 -l)C 2rfCo ' 



(49) 



and the integrals over 9 can be changed into integrals over ( . The boundary 
between long and short wavelengths is given by Co — 1 which corresponds to 
9 = 9 cr , given by [17]: 



sin 6 r 



M 2 



7 + 1 



2jM? + (3 - ^Mf - 2 



(50) 



We can decompose K as the following sum: K = 
using ( 149]) and taking into account Eqs.(l38]). ( 139 



K\,D + ^2D 



+ K$n. where 



42) and (H3 



"2D- 

we write: 



K. 



2D 



2 

7T JO 



9 roo 

\v rot \ 2 d0 = - 

IX Jl 



1 + 



:i-miko 2 

Ml 
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s 

rot I 



Mi 



i?(M x 2 



M?R + [Ml - 1)C 2 



dC 



o • 
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K. 



2 rl 2 



2D 



K ac 

-"-2D 



7r Je, 



v rot \ 2 d9 



7T JO 



;i - M|)c 2 



7T JO 



\v ac \ 2 d6 = - 



7T Jl 



Cr 



-1/ |2 
rot I 



M]_Jr{mi 



Mi J R(M 2 



M 2 R + (Mf 



ixf Co 



M 1 2 i?+(M 1 2 -l)C, 



2 ^Co 



In FigJHJa we plot the turbulent kinetic energy generated downstream sepa- 
rated in its different contributions as a function of the shock Mach number. 
We observe that the limiting value for the strong-shock limit is a function 
of 7. In FiglSJb we show the total kinetic energy generation as a function of 
7 and Mi, in the ranges 1 < 7 < 2, 1 < Mi < 10. It is observed that the 
total energy K 2 d grows unbounded in the strong shock limit (Mi — > 00) of 
a highly compressible gas (7 — > 1). 
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Figure 3: (a) Turbulent kinetic energy generated downstream separated in 
its different contributions as a function of the shock Mach number for an 
ideal gas with 7 = 5/3. (b) Total kinetic energy generation as a function of 
7 and Mi. 



(51) 



3.2 Acoustic energy flux 

The emission of sound by the shock is different in the long and short wave- 
length regimes, as discussed in Section 2 and with more detail in [T6| ITT] . In 
fact, for the long-wavelength regime, the acoustic field decays exponentially 
following the shock and its asymptotic contribution vanishes. However, the 
short wavelength modes contribute with stable running fronts downstream. 
We study here the acoustic energy flux in two different reference systems: the 
compressed fluid and the shock reference frames. The energy flux is defined 
by 03]: 

q = c 2 Ek ac , (52) 
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where E is the energy density of the sound wave: 

E = p 2 c\f , (53) 
in which p is the pressure fluctuation behind the shock and can be retrieved 



with the aid of Eq. ff24|) . The unit vector k ac can be expressed with the aid 
of Eq. (l27p . and represents the direction of the sound waves emitted by the 
shock downstream. We have: 

? k a x °x + k y y 

k ac = - == = (COSfln^Smfln 




2 



- M 2 

, ■ (54) 

,Co-M 2V /Co 2 -l Co-M 2V /Co 2 -l/ 

The projection of Eq. (l52j) along the longitudinal direction gives: 

q x = p 2 c\p 2 cos 9 ac = p 2 c\e\ cos ac . (55) 
For simplicity, we define the dimensionless longitudinal energy flux as: 

q x = — = e 2 s cos 6 ac , (56) 

where Co can be expressed as a function of the acoustic angle 9 ac 
r - 1 ( l-M 2 cos9 ac \ 

Co " { sine ac ) ■ (57) 



L 2 

In the new variable 8 ac , the region which corresponds to sonic waves traveling 



to the right is cos" 1 M 2 < 9 ac < vr/2 (1 < Co < l/Jl-Mi), and the 



corresponding to left-facing waves is ir/2 < ac < n (1/yl — M| < Co < °°)- 
We note that q x = for 6 ac = tt/2 which is trivial, and also vanishes for 



Co = Co = Mi/ J Mi — 1, because e s (Co) = 0. As has been done for the 
kinetic energy, we make the corresponding averages of the sonic flux over 
the dimensionless frequency Co in the compressed fluid frame. The acoustic 
energy flux can be separated depending on the direction of the traveling 
fronts: 



. ht rVy/^l M 1V / J R(M 1 2 - 1) 

{q * 9 )e = Ji qxiCo) MlR+{Ml-l)Q dCo ' 



f oo MiJR(Mf-l) 
= L<V^l UC0) M ? / + (M?-1)V C ° ' (58) 
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It is also interesting to express the sound energy flux in the shock reference 
frame, because in some experiments the shock remains steady with respect 
to the laboratory walls, and it is the fluid upstream that moves toward the 
shock. In that case, the normal to the sound wave front (which coincides 
with k ac ) does not coincide with the direction of propagation of the energy. 
The energy flux is now given by: 

q = v ac E s , (59) 

where v ac = c 2 k ac + (U — D)x, and E s is the energy density in the shock 
frame, it is written as H31: 



E s = (M 2 - cos 9 ac )E , (60) 
where E is given by Eq. ()53p . Collecting these results we obtain: 

q x = -^ = e 2 s (l - M 2 cos6 ac )(M 2 - cos6 ac ) . (61) 
P2C 2 

In the shock reference frame, the angle between v ac and the £-axis is denoted 
by 6'. Its relation to 9 ac can be seen to be given by: 

a> cos 6 ac -M 2 { 

cos 9 = —j^^^^^=^^^^^= . (62) 

x /l + M 2 2 -2M 2 COS Oac 
The averaged acoustic energy flux is given by: 



L ac 

IX JO 



q x d9 , (63) 



ikin = Pi^r = ^p2c\MlR i , (64) 



and we define the kinetic energy incident through the shock as: 

D*_ _ 1 

T ~ 2' 

The dimensionless emission coefficient S, which connects the incident kinetic 
energy flux with the acoustic energy flux emitted by the shock is: 

S = ^xel, (65) 

■* kin 

which is not difficult to particularize for the 2D case: 



2 2/°° M^R{Ml-l) 
S2D = nMi&J 1 ^^MFl^' (66) 
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In FigJUa we plot the relative acoustic energy flux S 2 d as a function of the 
shock strength for an ideal gas with 7 = 5/3. In the strong shock limit 
(Mi — > oo), S2D reaches an asymptotic value dependent on 7. In FigJUb we 
show 3D-plot for the same quantity S2D as a function of the shock strength 
Mi and the adiabatic exponent 7. 




Figure 4: (a) Relative acoustic energy flux S 2 d as a function of the shock 
strength Mi for an ideal gas with 7 = 5/3. (b) Relative acoustic energy flux 
S2D as a function of the shock strength Mi and the adiabatic index 7. 



3.3 Density amplification downstream 

An important point regards the amplification of the pre-shock density per- 
turbation field due to the shock compression. Similarly as with other pertur- 
bation quantities, the density perturbation field downstream can be decom- 
posed as the superposition of a steady and non steady contribution: 

p(k y x, t) = p ac (k y x, r) + p en (k y x) , (67) 

where the acoustic term is given by Eq. (1281) : 

pac(k y x, t) = p(k y x, r) , (68) 

and the entropic term can be obtained with the aid of Eq. (fT3|) after subtract- 
ing the acoustic contribution: 



WW ) P '\ M 2 J Pi ' 

We define the average of the ratio between the downstream over upstream 
asymptotic density field as follows: 

(5ft) 1 (6 pi) 



(5p{) & (5p\) 
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G . (70) 



As for the dimensionless function G, we particularize for the 2D case, and 
explicitly calculate the long/short wavelength and the acoustic contributions: 



G 
G 
G 



l 


2 


rir/2 




2D — 


TT 


J Ocr 


MfM$ 


s 


2 


fdcr 


"(1-MfMf) 2 


2D — 


7T 


Jo 


MfM$ 


ac 

2D — 


2 

TT 


r8 cr 

/ e 
Jo 


2 s d9 . 



2(1- 



2(1 -MlMl 
MlMl 
M\Ml 



MfMi 



e s + l 



eir + 1 



de 



(W 



(71) 



We plot these quantities for a gas with 7 = 5/3 in Figj5la, and we observe 
that G2D > G l 2D , in contrast with the behavior in Figj3j where K 2D < K l 2D . 
The total density amplification is shown in a 3D-plot in Figj5jb as a function 
of 7 and Mi . 




Figure 5: (a) Density amplification downstream separated in its different 
contributions as a function of the shock Mach number for an ideal gas with 
7 = 5/3. (b) Density amplification G2D as a function of the shock strength 
Mi and the adiabatic exponent 7. 



3.4 Vorticity generation downstream 

We focus now on the vorticity generated by the shock oscillations. This 
vorticity is conserved for each fluid element after shock compression in the 
absence of viscosity. Therefore, in the compressed fluid frame the vorticity 
is steady. Using (I3~TI) . we have for any single- mode perturbation: 

U z = = -=r = il 2 p s + iZ 3 , (72) 

kyC2 \k\c2 sin. Pi 
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where S7 2 and f2 3 have been defined in Eqs. (j32p and (1331) . In order to ob- 
tain the corresponding averages, it is convenient to use k = \k\ for the di- 
mensionless vorticity Q z . With this idea, we define this new vorticity as: 
&z = u z /(k c 2 ). Thus, this quantity is used to get the dimensionless fac- 
tor W, which quantifies the averaged vorticity generated by the corrugated 
shock: 

R 2 ) = W x e 2 , (73) 

Again, we particularize the averaged quantity W for the 2D case, and it can 
be also separated in its different wavelength regimes: 



W l 2D 



w, 



2 r/ 2 

7T 

2 



2D 



7T JO 
•/ 



+4; 



+ 2Vt 2 Vt z ei r + ^ 



sin 2 Ode 



2fi 2 ^ 3 e s + O 



sin 2 Ode 



(74) 



In FigJEla, we plot W 2D , W 2D and W 2 d = W 2D + W 2D . In FigElb, we show the 
total vorticity generated behind the shock, for different values of 7 and M\. 
We see that W 2 d is divergent for Mi 3> 1 and 7-I < 1. This is because of 
the size reduction of the eddies in that limit. This divergent behavior would 
be smoothed in a real gas, because, viscosity effects would become important 
for strong shocks traveling in highly compressible gases. 









^ " (a) 


—W %D 






Mi 



W 2D (7, Mi) 




Figure 6: (a) Vorticity generated downstream separated in its different con- 
tributions as a function of the shock Mach number for an ideal gas with 
7 = 5/3. (b) Vorticity generation W 2 d as a function of the shock strength 
Mi and the adiabatic index 7. 



4 Summary 

An analytical linear theory for the interaction of a planar shock wave with 
a 2D random density field has been presented. The entropy spectrum is as- 
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sumed to be isotropic. First, the simpler problem of the interaction with 
a 2D single-mode perturbation has been developed, in which, the acoustic 
waves and the vorticity/entropy perturbations generated downstream have 
been obtained. Thanks to the isotropy of the upstream perturbations, the 
important statistical averages downstream can be easily calculated. The tur- 
bulent kinetic energy and vorticity generation behind the shock are obtained 
as a function of the shock strength M 1; and the gas compressibility 7. The 
acoustic energy flux (noise) emitted by the shock front is studied in the com- 
pressed fluid and shock reference frames. The amplification of the density 
perturbation across the shock is also discussed. The perturbation fields gen- 
erated downstream consist of vorticity/entropy structures which induce mass 
and momentum redistribution and its associated mixing process. The results 
shown here can be used to deal with other types of interactions (random 
acoustic field or isotropic vorticity perturbations). 3D isotropic spectra can 
also be studied with the tools shown here, by changing the probability density 
function for the 3D case as done in [17] . The model presented here can also 
be extended to account different boundary conditions, such as the presence 
of a piston, a free surface or an ablation front driving the shock. Besides, it 
might also be applied to study the effect if a second shock traveling into the 
compressed spectrum generated by the first one. 
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